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COMPUTING  EQUILIBRIUM  PATHS  OF  CONTINUOUS  FAMILIES  OF  ECONOMIES: 
A HOMOTOPY  TECHNIQUE  AND  ALGORITHMS 

By 


This  study  is  in  two  parts.  The  first  part  develops  algorithms  to 

underlying  technique  — to  be  called  the  homotopy  technique  — appears 
to  be  of  Independent  interest. 

The  homotopy  technique  suggested  here  is  related  to  the  homotopy 
principle  of  Eaves  and  others.  In  the  context  of  studying  continuously 
changing,  parametric  systems  an  idea  similar  to  the  homotopy  principle 

0 s t £ 1)  a homotopy  path  <f>:  [0,1J  — * M “ [0,1)  is  traced  such  that 
$(t)  approximates  a solution  of  F(*,t)  for  all  t e [0,1).  Starting  at 

using  complementary  pivoting  until  it  terminates  with  a solution  to  the 
final  system  F(',l).  The  starting  solution  of  F(',0)  can  be  computed 


regarded  as  1 

the  homotopy  principle  with  a path  constraint  in  the 

following  ser 

ise.  The  path  initiated  at  the  trivial  solution  is  required 

to  join,  at  < 

in  intermediate  stage,  the  solution  path  of  the  given 

family  of  sys 

items  to  be  computed  and  follow  it  to  termination. 

equations,  g: 

ime  and  economic  equilibria  and  fixed  point  computations. 

But  in  this  v 

,«>  a,.  ...  1.  - equillbrlun 

version  of  Scarf's  algorithm  in  Che  sense  of  the  homotopy  technique. 


A new  kind  oi 

: nonsimplicial  subdivision  is  constructed;  its  cells  are 

' pseudocones 1 

in  g"  « [0,11  generated  by  ana  priori  chosen  subdivision 

of  S « [0,11 

.i  ...i.  1.  m,  .i  amaau  au  u •• 

subdivisions 

of  r”  x [0,11.  A refining  version  of  Scarf's  algorithm  is 

obtained  as  a 

i by-product. 

categories  of 

economies  considered,  a neoclassical  model  of  exchange. 

a general  Walrasian  model  with  production  and  more  elaborate  economies 
characterized  by  point-to-set  maps.  Each  of  the  algorithms  traces  an 
approximate  general  equilibrium  path  as  the  initial  economy  is 


be  adapted  fo 

r different  scenarios  like  (i)  only  the  initial  and  final 
specified,  (It)  the  initial  economy  and  a projected  path 

of  evolution 

are  given.  We  obtain  the  homotopy  principle  for  computing 

unknown  equilibria  of  a given  economy  by  taking  the  initial  economy 


implementation  i! 
subdivision.  Thi 


estimated  in  relation  to  fi  the  mesh  of 
approximate  path  converges  uniformly  t< 


Using  suitable  adaptations  of  the  main  algorithms  a new  dynamic 
framework  is  formulated  in  Part  II  for  the  analysis  of  various  problems 
in  economic  policy  evaluation,  public  goods  economics,  international 

1 prior  work  in  these  areas  using 


a has  been  done  in  a comparative  statics 
a dynamic  framework  had  been  emphasised. 


HOMOTOFY  TECHNIQUE  AND  ALGORITHMS 


CHAPTER 


INTRODUCTION  A 


' SUMMARY  OF  RESULTS 


This  study  Is  in  two  parts, 
of  formulating  algorithms  to  trat 
general  equilibrium  models  of  the 


i paths  of  approximate  equilibria  of 


The  method  underlying  these  algorithms  — to  be  called  the  homotopy 
technique  — which  provides  as  a specific  case  a parametric  version 
of  the  Scarf  algorithm  [107J  seems  to  be  of  Independent  interest. 

complementarity,  optimisation,  fixed  point  computations,  game  and 


general  equilibria  of 


algorithms  of  Part  I.  Problems  in  evaluation  < 
national  trade,  public  goods  economics  and  urbj 
new  dynamic  framework  for  analysis  using  paths 
the  underlying  systems. 

1.1.1  Homotopy  technique: 

A well-known  technique  for  studying  systems,  say  f:  M c R — f 

where  F(*,0)  is  a trivial  system  with  a known  solution  and  F(*,l)  ■ 
is  the  system  under  study.  Then  a homotopy  path  is  traced  starting 


n solution  and  terminating  in  a solution  of  F(*,l)  * f 


approximation  H of  F with  i 
a subdivision  of  M x [0,1] 
poiyhedra  such  that 
There  are  many  algorithms  whit 
45,  67,  91]  on  triangulations 


a subdivision  of  M*  (0,1]  Ci.e. , 
)-diraensional  closed  convex 
of  H to  each  polyhedron  is  affine), 
ov  such  paths  [26,  30-32,  34-35, 
(0,1]  by  complementary  pivoting. 


The  homotopy  principle  as  presented  by  Eaves  [32J  unifies  many  of  these 
algorithms.  The  emphasis  of  these  algorithms  and  the  homotopy  principle 
has  been  that  of  starting  from  a trivial  system  and  terminating  with  a 
solution  to  the  system  of  interest.  However  in  the  context  of 
studying  continuously  changing  (parametric)  systems,  essentially  a 
similar  idea  can  be  implemented  with  a different  emphasis. 

When  F:  M * [0,1]  — R represents  a one-parameter  family  of 
systems  (i.e.,  F(*,t)  for  0 s c £ 1)  a homotopy  path 


4>:  [0,1]  —Mx  [0,1]  can  be  traced 

with  a solution  to  the  final  system 

principle,  if  necessary. 

The  adaptation  we  have  made  wil 
view  of  applications  to  study  parametric 
discussion,  the  adapted  version  will  be  called  the  homotopy  technique. 
Though  at  first  glance  there  seems  to  be  no  essential  mathematical 


h that  $(t)  approximates  a 


he  following  sense, 
solution  path  of  the  given  family 


technique  ai 
principle  with  a patl 
starting  from  the  tr: 
intermediate  stage. 

Three  algorithms  are  developed  in  Chapter  5,  one  for  each  class 
of  general  equilibrium  models  discussed  in  Chapter  3.  The  first  two 
algorithms  compute  paths  of  e-approximate  Brouwer  and  Kakutani  fixed 
points  respectively.  The  third  algorithm  required  a distinctive 

ana  priori  chosen  subdivision  of  S x |0,1J 
advantage  of  this  type  of  subdivision  over 

subdivision  of  S « (0,1]  can  be  directly  related 


approximation  obtained.  As  a by-product  of  this  technique  a 
continuously  refining  version  of  Scarf's  algorithm  is  obtained. 

Chapter  3 describes  three  models  of  the  economy  of  varying 
generality  with  corresponding  characterizations  of  general  equilibria. 

point-to-set  maps  and  their  equilibria  are  characterized  by  versions  of 
the  Kakutani  fixed  point  theorem.  The  Walrasian  general  equilibrium 
model  with  production  is  treated  as  the  third  category.  Chapters  4 and 
5 develop  the  algorithms  for  computing  approximate  equilibrium  paths  of 
the  above-mentioned  classes  of  economies  under  continuous  deformation. 

The  homotopy  technique ‘is  presented  in  Chapter  4;  it  is  contrasted 
with  the  homotopy  principle.  The  underlying  tools  for  the  algorithms 

initial  economy  is  deformed  into  the  final  economy.  If  the  dynamics 
of  change  between  the  initial  and  final  economies  are  specified,  then 
the  algorithmic  output  approximates  this  temporal  evolution.  If  the 
dynamics  is  unknown  or  unspecified,  and  only  Che  initial  and  final 

paths  between  the  initial  and  final  economies.  If  the  Initial  economy 
an  approximate  path  of  evolution  of  the  equilibria.  As  a by-product 


technique  we  obtain  an  algorithm  to  compute  equilibria  of 


IV o kinds  of  results  obtained  formalise  the  sense  of  the 
approximation  implied  in  the  algorithms.  The  first  is  a limiting 
argument  which  asserts  that  as  the  mesh  of  the  simpliclal  subdivision 
is  reduced  to  sero  the  approximate  equilibrium  path  converges  Co  the 

In  Part  II  the  algorithms  of  Chapter  5 are  applied  to  study 
various  problems  In  economics  based  on  general  equilibrium  paths  of 
the  underlying  systems  subject  to  continuous  deformations  induced  by 
the  issue  under  study.  A new  dynamic  framework  has  been  devised  to 
study  dynamic  paths  of  policy  alternatives  and  economic  reforms. 

The  tax  structure  allows  dynamic  modeling  of  various  differential 

Chapter  7 develops  a similar  dynamic  framework  for  analyzing  issues 
in  international  trade.  The  framework  allows  for  a dynamic  analysis  of 


d by  changes 


international  markets  under  continuous 
markets  and  ocher  issues  in  world  trade, 
in  urban  economics.  The  effects  of  continuous  changes  in  urban  trans- 
area (4-5,  79-81)  has  used  a static  model.  Arnott  and  MacKinnon  (5, 

models  under  deformation. 

specifically  discussed  in  this  paper.  Prospective  applications  include 
parametric  versions  of  various  problems  so  far  tackled  using 
complementary  pivot  theory  (see  Section  1.4  and  Remark  1.3.4). 


1,2  The  Economic  Setting 

century  has  been  that  of  general  equilibrium.  Its  major  theme  is  that 
rium  arising  from  the  interaction  of  a number  of  economic  units  with 

production  of  a variety  of  consumer  goods  and  industrial  products, 


should  generate  a si 


es  of  political,  administrative  ai 


I institutional 


The  general  equilibrium  has  been  studied  at  different  levels  of 
generality  (6,  11,  17-25,  41,  54-58,  74,  107-120,  122,  127J;  the 

ness  and  stability  — all  these  issues  have  been  studied  [6,  8,  11,  17, 
18,  21,  23-25,  42,  109,  116-120,  122,  127],  and  continue  to  be  studied 


[97-99],  externalities  of  various  kinds  [11,  20, 
differences  in  agents'  information  [6,  58,  122], 

tions  under  which  an  equilibrium  exists,  sustains  its 
maintains  its  stability  and  so  on. 

The  main  underlying  devices  used  in  the  study  of 
systems  arising  in  general  equilibrium  models  have  be 
point  theorems  from  Algebraic  Topology.  Brouwer's  fi 

tions  and  extensions  have  'provided  proofs  of  existence  of  si 

lost  three  decades  [6-8,  11,  17-25,  54-57,  128].  But  it  was  only  in 
Che  last  decade  that  these  existence  results  of  general  equilibrium 

so-called  fixed  point  algorithms . The  pioneering  work  of  Scarf 


.5,  122,  130], 
the  equilibrium 

Its  optimality, 
of  the  complex 


[12. 


. 123-12i 


131-132]  has  made  it  possible  to  compute  the  equilibria  of  fairly 
large-scale  economies  with  reasonable  effort.  This  possibility  of 
obtaining  computed  values  for  general  equilibria  has  opened  up  a new 
framework  for  studying  numerous  problems  using  models  which  incorporate 
all  the  complex  interrelationships  of  general  equilibrium  theory  [4-5, 
46,  49,  54,  78-81,  97-99,  107-115].  Some  of  the  areas  studied  in  this 
framework  have  been: 

Analysis  of  taxation  policies  [5,  51-53,  108-113,  115], 
Impact  of  price  distortions  [114,  130], 

Problems  in  urban  economics  [4-5,  79-81], 

Public  goods  economies  [97-99], 

Issues  in  international  trade  [46,  54,  114,  129-130]  and 
A general  class  of  spatial  equilibrium  models  |80-81]. 

All  these  and  similar  other  studies  have  used  a comparative  statics 
procedure,  i.e.,  they  have  compared  Che  equilibria  in  the  situations 
before  and  after  the  change  under  the  study  was  made.  Clearly  such  a 
procedure  is  very  limited  in  that  it  provides  no  information  about  the 

etc.  The  thrust  of  this  paper  is  to  devise  algorithms  which  approximate 
equilibrium  paths  of  changing  economies  with  different  characterizations 
of  equilibria.  Some  applications  of  such  algorithms  to  areas  (1.2. 1-6) 
mentioned  earlier  are  also  explored. 


1.3  The  Problem  and  the  Approach 

The  research  in  this  paper  is  directed  at  achieving  three  major 


objectives: 


the  equilibrium  set  of  E 
n algorithm  I 


:) . The  problem  studied  here  is 
i connected  path  q(')  of  approximate 


equilibria  of  E(*).  the  sense  of  approximation  to  be  made  precise.  The 

principle  of  Eaves  [32]  and  others  with  a different  emphasis  in 
implementation.  This  technique  underlies  the  algorithms  of  1.3.1. 

lased  on  computed  solutions  in  a general 


In  the  last  eight 

the  algorithms  of  <1.3 
problems  by  providing  a dynamic  framework  ol 
opposed  to  the  existing  one  of  comparative  : 
The  first  two  objectives  are  discussed 


context  of  equilibria 


75-86,  89-96,  105-115,  131]  whose  solutions  have  been  attempted  using 
complementary  pivot  theory . 


Related  Approaches 

Our  development  of  the  homotopy  technique  and  the  algorithms 
Part  1 are  based  on  complementary  pivot  theory;  in  the  underlying 


ideas  and  its  presentation 
used  by  Eaves  in  (32J.  Hi: 
various  applications  of  a r 
complementary  pivot  theory 
The  applications  of  tf 
pervade  an  impressive 


t is  the  closest  to  the  homotopy  principle 
use  of  PL  Homotopies  in  [32]  unifies 
w and  poverful  technique  called 
hich  has  emerged  over  the  last  decade, 
e technique  in  its  various  forms 
f topics  including  linear  and  nonlinear 
t theory,  linear  and  nonlinear  optimization , 
equilibria  of  games  and  economies,  nonlinear  systems  of  equations, 
saddle-point  problems,  boundary  value  problems,  structural  mechanics, 
geometry,  plasticity  analysis,  network  problems,  pricing  of  utilities, 
urban  transportation  models  and  many  others. 

the  papers  of  Lemko  [75]  and  Lemke  and  llowson  [77]  where  a new 


underlying  principle  to  approximate  the  Brouwer  fixed  points  of  maps 
thereby  extending  the  purview  of  applications  of  this  combinatorial 


12 

technique  to  nonlinear  problems.  The  generality  of  the  underlying 
combinatorial  structure  was  recognized  [26-27,  67]  and  its  applications 

tial.  Linear  Complementarity  [13,  16,  28-29,  43,  62,  68,  70,  75-77, 
82-86,  95],  Nonlinear  Complementarity  [26,  37,  39,  60-61,  66,  89-90, 

Equations  (1-3,  10,  12,  27,  31-35,  40,  44,  48,  67,  71-73,  85,  92, 
122-125,  131-132]  have  been  probably  the  most  active  areas. 

Eaves  [30,  31],  Eaves  and  Saigal  [34]  used  the  idea  of  PL 
Homotopies  for  computing  fixed  points;  the  hemotopy  principle  was 

various  applications  of  complementary  pivot  theory  in  a crystallized 

this  dissertation. 

been  suggested  and  used.  Continuation  methods  [64-65,  96,  128],  Solving 
appropriate  differential  equation  systems  (65,  92],  Global  Newton 
methods  [120]  are  some  of  the  alternatives  which  can  be  used  if  the 
systems  under  study  satisfy  required  differentiability  conditions.  In 

are  made  on  the  characterizing  functions;  hence  we  use  simplical 


CHAPTER 


PRELIMINARIES  AND  SOME  USEFUL  THEOREMS 


Notation 


For  convenience  some  basic  definitions  and  conventions  for 


notation  and  referencing  are  collected  in  this  section. 


Enclidean  space  Rn  we  mean  Che  point  £ x^  where  £ Aj  - 1,  A^  e R; 

Given  a set  C then  aff(C)  (conv (CJ)  is  the  smallest  affine  (convex)  set 
containing  C.  Thus  aff(C)  is  the  affine  hull  of  C.  and  conv(C)  is  the 


2.1.2  Abbreviations: 

The  following  further  abbreviations 
bd  C ■ boundary  of  C. 


int  C ■ interior 


f C ■ dim  (aff(C) 
Matrices! 


14 


IF  A e R is  a real  m * n matrix,  we  denote  its  submatrices 

formed  by  the  rows  indexed  by  y in  their  natural  order.  Similarly  if 
B c {1,  2,  ....  n}  then  Ag  represents  the  submatrix  formed  by  the 


denotes  the  transpose  of  matrix  A. 

2.1.4  Set  Operations: 

Hie  usual  set  theoretic  operations  of  union,  intersection  and 
difference  will  be  denoted  by  u,  n and  ~ respectively. 


C + D - {x  + y e Rn;  x e C,  y « D) 
C - D ■ {x  - y e Rnj  x e C,  y e D). 


AC  - {Ax:  x e Ch 


C « D - {(x,y)  e Rn  + “j  x s C,  y s D). 


II  “ II 


B(x,p)  * Open  ball  with  center  x,  radius  p; 

l.e.,  {y  £ Rn:  ||  x - y ||  < p). 

B(C,p)  ■ p-neighborhood  of  C ■ C + pB(O.i). 
diam  C = Diameter  of  C » sup{||  x - y || : x,  y e C) 


15 


mesh  6 • mesh  of  G = sup(dl.lm  C:  C e 6) . 

where  6 is  a family  of  subsets  of  R . 


u is  said  to  be  lexicographically  lesB  than  v (denoted  as  u v) 


In  Chapter  3 the  preference  relation  for  i-th  consumer 
is  defined  on  R by  x i,  y iff  x is  at  least  as  much  preferred  as  y by 


2.1.7  Referencing: 


item  k in  Section  j of  Chapter  i. 
Section  i.j  refers  to  Section  j 


2.2.1  Definition: 


y t H(x,6)"n  C then  H(y)  c B(H(x),r). 

If  h:  C — * Rn  is  a function  then  H:  C --*■  P(Rn)  with  !l(x)  ■ (h(x)} 

iff  {h}  is  U.S.C. 


mappings 


following. 


2.2.2  Theorem: 
such  that  x^  — * x* , 


H(xq)  c B(H(x*),  e/2).  Thus  y*  c B(H(x*),  e)  for  a 
H(x*)  Is  compact,  y*e  H(x»). 

Two  additional  properties  are  stated  without  f 
p.  116-120)  or  Merrill  [91,  p.  17-19)  for  details. 
2.2.3.  Theorem: 

If  H is  a U.S.C.  point-to-set  mapping,  then  H( 
any  compact  subset  Dec.  [H(D)  ■ u H(x)]. 

2.2.6  Theorem: 

If  H is  U.S.C. , then  G:  C -+  P(Rn)  defined  by 


“ conv(H(x))  is 


When  we  consider  continuous  changes  over  time  on  a model  of  Che 
economy,  closely  related  are  questions  whether  one  map  can  be  continu- 
ously transformed  into  another.  The  elaboration  of  such  ideas  forms 
part  of  Homotopy  Theory  [88,  121). 


(0  I t s 1)  such  that  fp  * f and  f ^ c g.  Now  to  characterize  a 


continuous  family  of  maps  f : X — ► Y,  we  can  consider  a single  map 
F:  X * |0,1]  — > Y defined  by  the  rule  F(x,t)  ■*  ft(x)  (x  e X,  t e (0,1]). 
When  we  say  that  the  maps  f form  a continuous  family  of  maps  we  merely 


2.3.1  Definition: 


to  g)  if  there  exists  a continuous  map  F:X  * [0,1]  — Y such  that 
F(x,0)  - f(x)  and  F(x,l)  - g(x)  for  all  x e X.  The  map  F is  said 


be  a homotopy  and  we  write  f ~ g (f  is  homotopic  to  g). 


convex  subsets  of  R using  similar  arguments.  Since  the  maps  charac- 
terizing the  economic  models  often  are  those  into  R or  its  convex 
subset^  Theorem  2.3.2  guarantees  the  existence  of  a homotopy  to  deform 
the  initial  economy  into  the  final  one.  Such  a homotopy  can  be  used 
(when  the  dynamics  of  change  from  the  initial  economy  to  the  final  one 

The  trajectory  of  fixed  points  of  continuous  functions  under  liomotoplc 


Chapter 


Some  Underlying 


Some  fixed  point  theorems,  useful  in  characterizing  equilibria  of 
below.  These  theorems  are  stated  in  a manner  best  suited  to  the 

2.4.1  Theorem  (Brouwer): 

fixed  point,  I.e.,  x < S such  that  x - f(x). 

Equilibria  of  the  Walrasian  exchange  economies,  e.g.,  [6,  p.  26] 

dealing  with  fixed  point  algorithms  contain  a constructive  proof  [27, 

30.  39,  63-65,  71,  73,  91,  105,  107,  124,  132]. 

2.4.2  Theorem  (Kakutanl  [591): 

F:  S -4-  C(S)  is  a U.S.C.  point-to-set  mapping.  C(S)  is  the 
collection  of  nonempty,  closed,  convex  subsets  of  S.  Then  there  exists 

Equilibria  of  a general  Walrasian  model  with  production  [123,  p.  65) 
and  more  elaborate  models  of  the  economy  [6,  11,  17,  22,  46,  78,  97,  122) 
can  be  characterized  using  Kakutanl  fixed  points.  Constructive  proofs 


Of  direct  relevance  for  this  work  are  the  parametric  versions  of 
the  above  theorems.  A special  case  of  a theorem  of  Browder  [10,  p.  186] 


2.4.3  Theorem: 

Let  F:  S x [0,1) 

set  C In  S x [0,11  meeting  S « (0)  and  S * Cl)  such  that  F(x,t)  - x Co 
each  (x,t)  e C. 

Constructive  proofs  of  this  theorem  are  given  in  Caves  [26,  32) 
and  Freidenfelds  [39,  401-  Freidenfelds  [401  has  generalized  the  abov 
theorem  to  a parametric  version  of  Theorem  2.4.2  and  provided  a 
constructive  proof  using  complementary  pivot  theory  • 

2.4.4  Theorem: 

Let  F:  S * [0,1J  CCS)  be  a U.S.C.  point-to-set  map  and  CCS)  as 
in  2.4.2.  Then  there  exists  a closed  connected  set  C c s x [0,11 


« (0)  a 


« (1)  a 


technique  used  in  the  algorithms  of  Chapter  5. 


CHAPTER 


E ECONOMIC  FRAMEWORK 


Introduction  ' 

varying  levels  of  generality  are  discussed  below.  The  first  is  an 
abstract  model  of  an  exchange  economy  whose  equilibria  can  be  charac- 

more  elaborate  economies  where  the  defining  functions  are  point-to-set 
maps  and  their  equilibria  are  characterized  using  different  versions  of 

[107,  p.  76]  (to  be  referred  to  as  Scarf's  theorem  in  the  rest  of  this 
paper)  can  be  used  to  construct  the  equilibrium  configuration  of  a 
third  model  — the  Walrasian  general  equilibrium  model  with  production. 
In  Chapter  5,  three  algorithms  are  devised  to  trace  the  equilibrium 
paths  of  the  above  three  models  of  economies  under  deformation. 


In  this  section  we  consider  an  abstract  model 


good  or  service  completely  specified  physically. 


temporally  and  spatially.  It  is  assumed  that  they  are  perfectly 

The  commodity  space  can  be  represented  as  Rn,  by  stipulating  the 

i.e.,  inputs  (outputs)  are  represented  by  nonnegative  (nonpositive) 
real  numbers.  The  existence  of  all  futures  markets  will  also  be 


3.2.1  Preference  relation: 

Consumer  i has  an  initial  endowment  of  the  commodities,  w,  e R . 
He  also  has  a preference  relation  which  dictates  his  decision  when 
confronted  with  the  choice  between  two  commodity  bundles  say  x and  y; 
if  r y then  x is  at  lease  as  much  preferred  as  y.  Now  is  a 
binary  relation  defined  on  1C  which  is  reflexive,  transitive  and 
complete.  It  may  not  be  antisymmetric . Define  x y iff  x y and 


3.2.2  Strong  monotony  assumption: 


Bj(p)  = {x  e r":  pT(x  - Wj)  < 0).  If  p » 0,  Bj(p)  is  compact.  Clearly 


standard  simplex  S - (p:£p 
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If  the  preferences  further  satisfy  the  condition  that  for  all 
z r R,  , the  sets  lx  t R+:  z £,  x)  and  {x  e R+:  x £,  z)  are  closed, 
then  there  exists  a continuous  utility  function  u : r"  * R such  that 
u (x)  z ul(z)  iff  x z,  see  Debreu  117,  p.  54].  Such  preferences 
are  called  continuous  preferences.  When  is  continuous  then  the 


djfp)  - <x  c R":  x is 
Identity  p1d,(p)  ■ p for  all 


L<p)>- 

> 3.2.2)  we  have  the 


demand  function.  Then  p g(p)  * 0,  i.e.,  the  value  of  excess  demand 
is  zero  for  all  p c S.  This  equality,  called  the  Walras  law  is  a basic 
relationship  of  Walrasian  economics  [127].  Since  each  d,  is 
continuous,  g:  S — ► Rn  is  continuous.  This  exchange  economy  with  m 


agents  can  be  characterized  by  an  m-tuple  (~j,  w ) or  by  the  more 
derived  specification  (d^,' w() . On  an  aggregate  level  the  economy  can 


Id  g(?)  i 0. 


system  p 
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3.2.4  Definition: 

Note  that  If  Walras  law  holds  and  p Is  an  equilibrium  price  vector 
then  g (p)  ■ 0 for  all  l ■ 1,  2,  ....  n with  equality  if  p^  > 0.  Thus 
all  markets  are  in  balance  except  possibly  those  of  free  goods,  which 

(d( t j ) , w(tjO)  at  time  tj , the  algorithm  devised  in  3.2  will  trace 
the  equilibrium  path  of  the  economy  under  deformation. 


' More  Elaborate  Economies 


and  issues  of  welfare.  For  example,  see  Negishi  [96),  Debreu  [17, 
p.  37-881  or  Arrow  and  Hahn  [6,  p.  52-129).  In  these  models  the 
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economies  has  been  the  Kakutanl  fixed  point  theorem  and  its  variations 
applied  to  an  appropriate  U.S.C.  point-to-set  map  f:  S — * C(S),  where 


it  could  be  a simplex  of  utility  weights,  as  in  Neglshi  (96)  or  a 
simplex  of  normalized  vectors  of  utility  levels  of  all  the  agents  in 
the  economy,  as  in  Arrow  and  Hahn  (6)  or  some  similar  construction.  In 
some  cases  the  original  economy  characterized  by  maps  defined  on  an 
arbitrary  compact  set  C might  be  shifted  onto  the  standard  simplex  S 

denote  the  argument  of  f:  S ► C(S)  by  x instead  of  p. 

These  economies  are  so  varied  and  numerous  that  they  will  not  be 
described  in  any  further  detail.  For  our  purposes  the  functions 
f:  S — ► C(S)  shall  characterize  these  economies.  In  this  framework, 

fixed  points  of  the  maps  f(',t):  S — ► C(S)  as  t goes  from  its  initial 
algorithm. 


4 Production 


Now  we  discuss  the  general  Walrasian  model  obtained  by  adding  an 
explicit  production  sector  to  the  exchange  model  of  3.2.  As  before 
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each  consumer  will  determine  his  consumption  plan  by  maximizing  utility, 
subject  to  the  constraint  that  expenditure  at  the  prevailing  set  of 

inputs  and  outputs  being  evaluated  at  the  prevailing  prices.  The 
price  system  in  the  economy  forms  Che  link  between  these  independent 

that  no  producer  has  a compelling  motivation  to  search  for  higher 

As  in  3.2  w.  represents  Che  vector  of  resources  of  the  i-th 
consumer  and  d,(p)  his  demand  at  prices  p;  i.e.,  d.(')  is  the  demand 

analysis  matrix 


and  required  as  input  (if  Aj  < 0). 


activity;  if 


activity  j is  used  at 
as  output  (if  A^  2 0) 
first  n columns 


correspond  Co  an  assumption  of  free  disposal.  Hie  activities  can  be 

that  the  set  of  production  plans  available  to  the  economy  is  given  by 
the  set  {Ay:  y ■ 0}.  The  model  can  be  generalized  to  allow  for  more 
general  production  sets,  if  the  customary  convexity  assumptions  are 
made;  see  for  example  [17„  p.  39].  We  make  the  following  assumption. 


of  activity  levels  that  giv< 
1 commodities  is  a bounded  si 


Aggregating  the  individual  demand  functions  d <•)  we  r 
the  market  demand  function;  similarly,  w the  vector  of  tot. 

on  S,  as  before.  Ualrae  law  also  holds;  pTd(p)  = p w at  a] 
3. A. 2 Definition: 


e to  a nonnegative  net 
.e.,  {y:  y i 0,  Ay+w  l 0) 


(a)  supply  equals  d 


From  (a),  p w + p Ay  = p w,  and  Walras  law,  we  have  *p  Ay  ■ 0. 

Then  (b)  implies  p A.y^  ■ 0 for  j ■ 1,  2 k.  pV  = 0 when 

yJ  > 0.  The  intreprctation  is  that  no  activity  makes  a positive  profit 
at  the  prices  ; those  activities  used  (at  a positive  level)  must  make 

used,  the  price  of  the  corresponding  commodity  is  zero.  In  this  sense, 
at  the  equilibrium,  production  is  consistent  with  profit  maximization. 


be  based  on  versions  of  Kakutanl  fixed  point  theorem  (see  Debreu  (17, 
p.  831,  Scarf  (107,  p.  1201  or  Todd  [123,  p.  651).  Scarf's  algorithm 
of  [105  Or  107]  based  on  Scarf's  theorem  provided  the  first  construc- 


by-product. 

function  of  the  parameter  t,  for  tjSIS  tj,  where  (d( • , t) , w(t),  A(t)) 
characterizes  the  economy  at  time  t.  The  third  algorithm  we  have 
devised,  in  5.6,  will  trace  an  approximate  equilibrium  path, 

(p,y>:  [t0,t1l  -h-  S X R*  such  that  (p(c),y(t))  approximates  the 


The  model  of  this  section  provides  the  best  framework  for  the 


applications  of  Chapters  6 and 


BUILDING  BLOCKS  0: 


4.1  Homotopy  Technique 

The  underlying  machinery  of  the  algorithms  — which  we  call  the 
homotopy  technique  is  as  follows.  The  changing  system  under  study 
is  formulated  as  a one-parameter  family  of  systems  — a homotopy  which 
represents  the  actual  deformation,  and  at  the  initial  and  final  levels 

level  by  the  method  of  complementary  pivoting.  In  this  sense  the 

[32].  In  the  special  case  when  the  initial  system  considered  is  a 
convenient  trivial  system  homotopic  to  the  actual  system  under  study 
then  the  homotopy  technique  reduces  to  the  homotopy  principle. 

points,  zeros,  etc.)  is  to  imbed  f into  a one-parameter  family  of 
systems,  a homotopy  F:  H * [0,1]  — * Rn  where  F(*,0)  is  a trivial 
system  with  a known  solution  and  F(*,l)  = f is  the  system  of  interest. 


Then  a homotopy  path  is  traced  from  the  known  solution  terminating  in  a 
solution  of  F(*,l)  = f-  When  F does  not  satisfy  differentiability 


conditions,  the  strategy  of  tracing  the  homotopy  path  is  to  track 
solutions  to  piecewise  affine  approximation  H of  F with  respect  to  a 
subdivision  of  M * 10,1]  (i.e.,  a subdivision  of  H x (0,11  into 
(n+l)-dimensional  closed  convex  polyhedra  such  thac  the  restriction  of 
H to  each  polyhedron  is  affine) . There  are  many  algorithms  which 
follow  such  paths  on  triangulations  of  M * (0,1]  by  complementary 
pivoting.  The  homotopy  principle  as  presented  by  Eaves  [32]  unifies 
many  of  these  algorithms.  The  emphasis  of  the  homotopy  principle  has 
been  thac  of  starting  from  a trivial  system  and  terminating  with  a 
solution  to  the  system  of  Interest.  However  in  the  context  of  studying 
continuously  changing  systems  more  or  less  the  same  principle  can  be 

systems  (F(*,t)  for  t c (0,1J)  a homotopy  path  $:[0,1]  —*  M x [0,11  can 
be  traced  such  that  <|»{t)  approximates  a solution  of  F ( ■ , t ) for 
t e |0,1J.  The  starting  solution  can  be  computed  using  a first  phase 
application  of  the  homotopy  principle. 


The  adaptation  we  have  made  of  the  homotopy  principle  to  formulate 
the  homotopy  technique  will  be  very  important  from  the  point  of  view  of 
studying  parametric  systems.  Though  at  first  glance  there  seems  to  be 
no  essential  difference  between  the  homotopy  principle  and  the  homotopy 
technique  suggested  here  there  is  one  crucial  difference.  The  homotopy 

principle  with  a path  constraint  in  the  following  sense.  The  path 
starting  from  a trivial  solution  is  required  to  join,  at  an  intermediate 


stage,  a solution  path  of  the  given  system  (which  is  to  be  computed)  and 

follow  (32)  for  a presentation  of  the  underlying  tools. 

The  homotopy  technique  can  be  thought  of  as  a generalization  of 

to  study  systems  under  deformations  and  to  solve  parametria  versions 
lich  complementary  pivoting  has  been  applied 
le  work  of  Eaves  (30-32),  especially  (32),  has 


Chapter  3. 


economies,  i.e.,  the  systems  t*  study  are  various  models  of  the  econor 
and  the  solutions  of  Interest  are  their  general  equilibria.  Conse- 
quently the  algorithms  of  the  next  chapter  are  developed  in  the 
specific  coni 
appropriate  I 
the  detailed 
framework  wot 
broader  contt 


algorithms  we  follow  the  presentation  of  (321  very  closely. 


4.2  Cells  and  Manifolds 

convex  polyhedra  as  components  of  subdivision  of  the  sets  we  work  with. 


In  this  context,  after  Eaves  132],  we  will  call  then  cells.  Although 
the  basic  properties  are  well  known,  we  include  them  here  for  fixing 

By  an  n-cell  we  mean  a cell  of  dimension  n.  If  an  n-cell  is  the 
convex  hull  of  (n  + 1)  points  we  call  it  an  n-sirrplex. 

The  interior  int  o and  boundary  bd  0 of  a cell  0 are  the  interior 
and  boundary  relative  to  aff(O). 

Let  t be  a subset  of  an  n-cell  0.  If  x,  y e o,  0 < u < 1,  and 
(1  - p)x  + by  € T imply  that*  x c t and  y e T then  T is  called  a face  of 
a cell  0.  Faces  thac  are  (n  - J)-cells  are  called  facets  of  the  cell. 

The  intersection  of  a line  {x  + liy:  u e R]  and  a cell  o is  called 
a chord  of  o.  A nonempty  chord  may  be  of  dimension  0 or  1.  By  a ray 


4.2.1  Definition: 

Let  G ^ it  be  a finite  or  countably  infinite  collection  of  n-cells 
in  Rn.  Let  G*  be  the  set  of  i-faces  of  the  elements  of  G.  G ■ u a. 

We  call,  (G,G)  a subdivided  n-manifold  if 


s a neighborhood  meeting  only  a finite  number 
re  is  a G,  such  that  (G.G)  is  a 


subdivided 


e subdivision  G ii 


4-2.2  Definition! 

For  (G,G)  os  in  4.2.1  if  each  n- 
n-simplex  we  say  that  G triangulateo  G. 

Of  particular  importance  to  us  in  constructing  equilibrium  paths 

curve.  It  is  trivial  to  prove: 

4.2.3  Lemma: 

A connected  1-manifold  (i.e.,  a curve)  is  homeomorphic  to  either 
a circle  or  an  interval.  For  the  two  cases  above  we  call  the  curve  a 
loop  or  a route  respectively. 


1-manifold  is  a disjoint  collection  o: 


Definition: 

be  an  n-manlfold  subdivided  by  G.  By  the  boundary  of  G, 
ean  the  union  of  all  (n  - l)-cells  of  G which  lie  in  exactly 


s (32,  p.  85]. 


4.2.8  Definition: 

Let  P be  a 1-manifold  neat  in  G and  let  6 be  a subdivision  of  G. 
n-cell  of  G.  If  P subdivides  P we  say  that  P is  neat  in  (G,G)  or 


4.3  Piecewise  Affine  Maps 
4.3.1  Definition: 


Let  (G,G)  and  (J,J)  be  subdivided  manifolds . Let  F : G — * _ 
a continuous  map  which  is  affine,  (F(px  + (1  - u)y  “ pF(x)  + (1 

Such  a map  F is  called  piecestise  affine  map. 


Given  a cell  a E G and  T e J with  F(a)  c t define 
F : aff(a)  — ♦ aff(T)  to  be  the  affine  map  which  agrees  with  F on  O. 

all  y c aff(G).  This  representation  is  assumed  to  be  with  respect  to 
the  standard  basis  of  R. 


4.3.2  Definition: 

affine  map  K:  C — * Rn  by  taking  K(x)  »Iu  H(x  ),  E p = 1,  p.  J 0 
called  a piecewise  affine  extension  of  H.  When  there  is  no  ambiguity 


manifold  in  Rn  let  o € G have  verticies  x_,  x.f  x^,  ....  x , (thought 

A = (H(x.)  - H(x_),  ...»  H(x  ) - H(X-))  be  n x n matrices.  Then 

Proof:  Let  x ■ x^  + Wp  e a where  p ■ (p*.  p , ....  pn) . Hence 
p - H_1(x  - XQ).  H0(x)  - H(x0)  + AP  - H(xq)  + A W_1(x  - X,,)  . 

4.3.4  Definition: 

which  agrees  with  H on  the  vertices  of  G (a  triangulation  of  G)  is 
called  a pieceuise  affine  approximation  to  H with  respect  to  the 
triangulation  G.  Given  such  an  H and  H,  we  hove: 

If  ||  x - x' ||  S ( implies  ||  H(x)  - H(x')|l  it  and  if  mesh  of 
G ; 6 then  ||  H(x)  - fi(x)||  i e for  all  x e G. 

In  following  the  paths  of  equilibria  using  complementary 


e algorithm  of  Eaves  [32]  is  ; 


e algorithms  of  the  next  chapter 
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boundary  start.  This  algorithm  forms  a b. 
algorithms  in  their  implementation. 

The  algorithm  is  based  on  1-manifolds  neat  in  (G,G)  which  are 
inverse  images  of  points  y € J under  a piecewise  affine  map  F : G — *• . 

F_1(y)  is  a 1-manifold  the  point  y has  to  satisfy  a weak  regularity 

f these  1-manifolds  in  (G,G) . 

Definition: 


Let  G,  J and  F be  as  in  the  paragraph  above.  A point  x in  G is 
said  to  be  a degenerate  (otherwise  regular ) point  of  F if  x lies  in  a 
cell  0 s G with  dim  F(o)  < n.  A value  y in  F(G)  is  defined  to  be 
degenerate  (otherwise  regular ) value  if  F (y)  contains  any  degenerate 


points  is  regular.  Clearly  for  a 


Eaves  algorithm. 

4.4.2  Theorem: 

If  y is  a regular  value,  then  F 1(y)  is  a 1-manifold  neat  in  (G  ,G)  , 
and  F ^(v)  is  subdivided  by  sets  of  the  form  0 n F * ( y ) it,  0 € G . 


regular  and 


boundary  points  plus 


(G,G)  and  <J,J)  a 
n respectively.  F:  G ■ 


(y)  1 


algorithm  of  [321  describes  a strategy 
')  containing  in  the  direction  u_. 


for  moving  along  a curve 
The  algorithm  termination  depends  on  Che  nature  of  the  curve  of  F (y) 
in  the  direction  u and  G.  The  algorithm  may  terminate  in  a finite 
number  of  steps  with  one  of  the  three  possibilities,  i.e.,  a "Ray," 


then  the  algorithm  does  terminate  in  a finite  number  of  steps.  Then  if 


d G or  a "ray"  is  guaranteed, 
important  part  in  the  algorithms  of  Chapter  5. 


Now  we  sketch  briefly  the  algorithm  ol 
s implemented  with  a boundary  start. 


r regular  values 


s a regular  value, 


L (as  in  4.3.3)  provide  a direction  uQ  ! 

g > o).  Thus  we  have  the  initial  triple. 
Iteration  r = 0,1,2,...:  From  the  triple 
6r  - sup{6:  x + 6u  e a ) and  record  (xf 
register  a "ray"  termination  and  stc 
If  xr+1  £ bd  G,  register  a boundary 
determine  e G ~ tor)  which  coni 


1 sufficiently  smai: 
' to  iteration  r » 0 


(r+1)  with  the 


Given  a G and  x a facet  of  an  (n 


0 £ G one  should  havi 
£ G - (o)  which  contj 
s,  called  pivoting 


commonly  used  triangulations. 


nature  of  G and  F.  In  the  algorithms  of  Chapter  5 steps  similar  to 
feasible  pivot  as  in  the  Simplex  method  of  Linear  Progratmning  (15). 


pivoting  ruins  of  familiar  triangulations  would 


a simplicial  subdivision. 

Note  that  Algorithm  4.4.6  wi 

is  specified  for  regular  values.  An 

elaborate  version  is  described  ii 

» [32]  for  handling  degenerate  values. 

Conceptually  the  idea  for  treating  a degenerate  value  y is  to  perturb 


it  to  a regular  value  y + Y[e]  wl 

.ere  fc ] is  a vector  (e,  ez,  ez e”) 

of  powers  of  e and  Y - [y1>  y2. 

. ...  y ] is  a matrix  of  rank  n.  At  the 

implementation  level  the  method  : 

is  to  use  the  lexicographic  technique 

for  resolving  degeneracies.  For 

details,  see  [32,  p.  111].  Similar 

algorithms  of  the  following  chapter.  Majthay  (82,  84]  and  Caves  [28] 


pivoting. 


CHAPTER 


THE  ALGORITHMS 

economy  is  changed  from  Its  initial  state  at  time  t-.  to  its  final 
state  at  t,  according  to  a specified  dynamics.  If  the  dynamics  are 
not  known  or  not  specified  then  any  homotopy  (see  2.3)  between  the 
two  economies  provides  possible  deformations  of  the  initial  economy 
into  the  final  one;  in  this  case  the  algorithms  trace  possible 

when  the  initial  economy  and  the  dynamics  of  evolution  is  specified 
but  the  terminal  point  is  left  free.  Then  the  algorithmic  output 
provides  a path  of  equilibria  on  the  semi-open  interval  [t_,t.)  where 


can  be  loosely  described  as  follows.  The  economy  w 


first  two  algorithms  and  K ' 
algorithm,  C(S)  (as  in  2.4.2 
third.  The  deformation  of  t 


algorithms 


F:  H k [tg.tjl  — * N.  The  function 

scarcing  at  Xq . The  path  will  terminate  in  M * {t^}  at  an  equilibrium 
of  the  final  economy.  This  is  guaranteed  by  remark  4.4.5  and  the 
formulation  of  the  problem  such  that: 

(a)  H L(y)  contains  no  "rays,"  and 

(b)  H_1(y>  n((bd  M)  * [tp.tj]}  - 4. 

The  above  strategy  underlies  each  one  of  the  algorithms  to  be 

ie  algorithms  and  will  depend  on  the  description  of  the  economy 
he  equilibria. 


description  of  the  algorithms  we  shall  always  consider  the  intervals 
of  a homeomorphism  the  functions  involved  can  be  reparametrized. 


Che  sense  of  die  approximation  obtained.  The  first,  foremost  in 
computational  considerations,  is  the  closeness  of  the  approximation 
of  the  equilibria  in  relation  to  tl 
second,  of  conceptual  importance,  : 
converges  to  the  actual  equilibria 

123-124,  131 ) . For  a comprehensive  discussion  of  triangulations. 


limiting  sense,  i.e.,  when  the 


complementary  pivoting. 


then  the  triangulation  will  be  called  a uniform  triangulation.  If 
fi.  -*  0 as  k — *■'»  it  will  be  called  a refining  triangulation.  Merrill 
[91J  Kuhn  and  MacKinnon  [731,  Kuhn  [71],  Gould  and  Tolle  [47]  and  Todd 
[1251  contain  versions  of  uniform  triangulations.  Eaves  130) , Eaves 
and  Saigal  [34]  and  Todd  [123-1251  contain  detailed  discussions  of 
refining  triangulations.  These  references  also  contain  specific  rules 


5.2  Algorithm  for  the  Exchange  Economy 

In  this  section  we  devise  an  algorithm  to  trace  an  equilibrium 


Recall  that  these  economies  could  be  specified  in  an  aggregate  form 
using  the  excess  demand  functions  g:  S Rn  and  their  equilibria 
could  be  characterized  (by  Definition  3.2.4)  as  p e S such  that 

g(p>  i o. 

representing  the  change  in  economies  from  g(*,0)  to  g(',l)  is  specified. 
(p,t)  would  be  an  equilibrium  at  time  t if  g(p,t)  £ 0.  We  construct  a 
map  f:  S x [0,11  — *■  S that  the  fixed  points  of  f(*,t)  coincide  with 

5.2.1  Construction: 

f:  S x [0,1]  — ► S is  defined  as  follows.  For  (p,t)  e S x (0,1) 
set  for  i - 1,  2 n 

f^p.t)  - J^-ax|0,B1<P,t>l 
1 + E max[0,g  (p,t)] 


where  f*  and  g*  denote  the  i-th  component  of  f and  g respectively. 

A fixed  point  of  f(',t)  for  any  t e (0,11  represents  the 

cp1  = p1  + maxIO^Cp.t)]  where  c = 1 + S max[0,g^(p,t) 1 . If  c > 1, 

£ 

Chen  pi(c-l)  - max|0,g1(p,t)l  Implying  that  g*(p,t)  > 0 for  all  i with 
Thus  c ■ 1,  and  g(p,t)  i 0. 


5.2.2  Remark: 

Thus,  by  construction  above,  the  problem  we  have  reduces  to 

given  an  e > 0,  the  algorithm  will  compute  a continuous  path, 

II  fWt»  - ♦jCDlU  E for  any  t e [0,1].  The  strategy  in  the 
algorithm  is  to  compute  a path  of  the  fixed  points  of  g:  S * (0,1J  — ► S 


||f(x,t>  - f(x',t)  ||  S e whenever  ||  x-x'  ||  £ 4. 


f(-,t)  for  all  t « [0,1], 


Proof : If  O(t)  is  a fixed  point  of  g(-,t)  we  have  g(a(t) ,t)  ■= o(t) . 

have  II  f(a(t),t)  - §(a(t),t)||  S 0.  Thus  ||  f(c.(t),t)  - a(t)  ||  i € as 
required. 


r f' 

p.o 


where  ( p,  t)  is  a vertex  of  G®  in  S x [0,11,  < p_ , 0 ) is  the  equilibrium 

S x (0}  is  artificially  specified  such  that  (p_,0)  is  the  unique 
solution  to  H(po,0)  " 0.  If  p_  is,  in  fact,  the  unique  equilibrium 

piecewise  affine  map  H:  S x [0,1]  — x S - S as  in  4.3.2. 

from  (Pq>0)  to  (Pj.l)  is  the  path  of  fixed  points  of  g,  and  thus  a 
path  of  s-approximate  fixed  points  of  f.  Thus  we  obtain  the  required 
path  of  approximate  equilibria.  The  details  of  constructing  this  route 


Remark: 


an  imbed  S in  a larger  simples  S'  and  artiflcally 
obtain  the  condition.  The  impossibility  of  boundary  equilibria  has 
also  been  evoked  based  on  economic  arguments  as,  for  example,  in 
[21-23,  95).  An  assumption  of  this  sort  will  be  made  in  all  the 

5.2.7  Remark: 

Provisions  for  handling  degenerate  values  are  incorporated  in  ou 
algorithms.  The  underlying  concept  is  to  perturb  the  starting  point 
(which  is  made  a regular  point  by  construction  of  H)  to  qn  such  that 

(lexicographic)  linear  systems.  In  details,  we  follow  Eaves  [32). 
Start:  Given  a triangulation  0 of  S * [0,1)  construct  H as  in 

0-  be  the  unique  (frt*l)-celL  containing  Y«.  To  ensure  the  regularity 


points  into  o0  from  qQ  + Qq ( E ] . (qQ,  (?0,  Op,  u0)  is 
quadruplet.  Note  that  [el  plays  no  real  part  in  the 
pivoting  in  the  corresponding  lexicographic  system  i 
H(qQ  + Q0 [e] ) can  be  shown  to  be  a regular  value  for 
With  (q„,  Qn»  On,  un)  as  input,  go  to  iteration  0. 


tb:  qr  + Q [el  + b u e or)  for  all  small 

[15)  or  [32,  p.  1101).  bet 
Vfl  + - qr  + (^[el  + (br  + Br[el)ur 

11  small  e > 0.  If  this  point  is  in  S « {1}  for  small  e > 0, 
Otherwise  determine  the  unique  r G - (o^l  which  contains 

points  into  Proceed  to  iteration  (r+1)  with  the  quadruplet 

* [0,11  is  not  even  considered  in  the  algorithm. 

number  of  iterations  when  G is  a uniform  triangulation. 

This  concludes  the  discussion  of  the  basic  algorithm.  Now  some  of 
its  properties  and  adaptations  for  solving  related  problems  are 


discussed 


If  the  economics  arc  specified  only  at  Che  initial  and  final 
levels,  and  Che  deformation  of  one  to  the  other  is  unknown  or 


unspecified  we  can  use  this  algorithm  to  compute  possible  equilibrium 
paths  connecting  the  equilibria  of  the  initial  economy  to  those  of 
the  former.  Assume  that  the  initial  and  the  final  economies  are 
characterised,  as  in  5.2.1,  by  f(’,0)  and  f(’,l):  S — * S.  Now  using 


K(-,0)  - f(-,0)  and 


The  algorithm, 
path  (of  equilibria) 


Implemented  using  some  such  F,  produces  a possible 
of  smooch  transition  from  an  equilibrium  of  the 


If  the  initial  economy  is  specified  along  with  f:  S > 10, “)  -•  S 

mesh  of  the  triangulaclon'of  S x (0}  is  taken  Co  be  6 chosen  as  in 
5.2.2.  Eaves'  trlnngulation  K(  and  K-  of  (30)  could  form  the  basis 

triangulation  of  5 ■ (0,1). of  mesh  6 (as  in  5.2.2).  This  would  Involve 
a straight-forward  reparametrization  x:  [0,™)  — » (0,1J. 
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5.2.12  Algorithm: 

The  algorithm  can  be  initiated  at  (p(0) , x(0))  * ( pQ , 1 > in  H-*(0) . 
The  algorithm  generates  a route  which  we  can  parametrize  by  (p(t),T(t)) 
where  t c [0.+  -)  and  (p«),T<£)>  - q£  c <J£  for  l - 0.  1,  2.  ... 

the  details  of  implementation  are  as  in  5.2.8.  Note  that  even  though 
the  algorithm  generates  an  infinite  sequence  of  simplices,  it  can  be 
terminated  at  any  level  t,  with  an  equilibrium  path 
*:  10, tj]  S x [O.tJ. 

This  alogrithm  can  be  used  to  compute  an  unknown  equilibrium 
point  of  a given  static  economy  f:  S — * S,  as  is  well  known.  Using  a 
refining  triangulation  of  S x (0,1]  and  the  homotopy 
f(p»c)  • t(p  - p_)  + (1  - t)  f(p)  where  (p_,l)  is  an  arbitrarily  chosen 
starting  point.  Algorithm  5.2.12  can  be  used  to  generate  a path 
starting  at  (p.. , 1)  to  an  equilibrium  point  of  the  given  economy  as 


to  S * to}  as  (p_,0)  is  a unique  solution  to  H(p,0)  ■ 0 


by  construction. 


the  path 


possibility  of  retrogression 


however,  this  retrogression  can  be  prevented.  The  strategy  is  as 
follows.  At  any  level  k,  k j*  0 let  0.  be  Che  n-slmplex  containing  the 
point  p,  such  that  H(p.  ,k)  ■ 0.  After  0.  has  been  encountered  for  the 
first  time,  use  at  the  k-level,  an  artificial  piecewise  affine  map  such 
that  (p.  ,k)  is  the  unique  solution  in  S * { k } (instead  of  the 
originally  defined  H) . A(p-p.)  provides  a convenient  such  map,  where 
A is  a suitable  nonsingular  matrix.  The  algorithmic  path  never  returns 


we  obtain  a connected  set  meeting  S * {0}  and  S ' (1).  providing  a 


characterized  by  point-to-set  maps.  The  underlying  strategy  and  many 


discussed  in  3.3.  In  these  models  the  functions  characterizing  the 
economics  are  point-to-set  maps;  their  equilibria  (as  discussed  in  3.3) 
can  be  taken  to  be  characterized  by  the  Kakutani  fixed  points  of  an 


appropriate  U.S.C.  point-to-set  map.  So  in  this  section,  for 


1 characterize  t 


points,  (x  t S:  x c f(x,t)}  coincide  with  the  equilibria  of  the  economy 
at  time  t.  In  these  models  S,  the  standard  simplex  in  an  appropriate 

economy  from  the  initial  one  f(',0)  to  the  final  one  f ( ’ , 1)  . Our 
algorithm  would  trace  an  approximate  equilibrium  path  starting  at  an 
equilibrium  (xfl,0)  of  the  initial  economy  and  terminating  at  an 
equilibrium  (x^ , 1)  of  the  final  economy.  The  underlying  strategy  is 
to  compute  a route  of  approximate  fixed  points  of  f.  More  precisely: 


Given  f:  S X [0,1)  — ► C{S)  such  that  f ( • , t)  characterizes  the 

rtj.'tj):  [0,1]  — S x [0,1)  such  that  *2<0)  = 0,  *2(1>  - 1 and 
en  any  E > 0,  ^(t)  « B(f(*(t),e)  for  all  t c [0,1|. 

Before  the  details  of  the  algorithm  are  discussed  we  need  the 


5.3.2  Definition: 

Let  f:  S — x PCS)  be  a point-to-set  map.  Given  e > 0,  x e S such 
that  x c B(f(x),e)  is  called  an  e-approximate  fixed  point  of  f. 

5.3.3  Definition: 

all  x € S.  For  each  vertex  y of  6,  pick  g(y)  e f(y)  and  let  f:  S — x S 
be  a piecewise  affine  map  which  assumes  the  preselected  value  g(y)  for 


defined 


called  a piecewise  affine  approximation  c< 
local  finiteness  of  G it  is  easy  to  show  i 


y s B(x,«)  n S,  f(y)  c B(f(x),€),  as  f is  U.S.C. 

Now  if  x is  a fixed  point  of  g,  then  x ”i|1  X^  Ajftyj) , 

g(yt)  « f(y4)  = B(f(x),E)  for  i = 1,  2,  ....  n.  Since  B(f(x),e)  is 
convex  it  contains  x ■ £ X.  y.,  fl,  ■ 1,  X,  2 0,  i.e.,  x e B(f(x),e) 
is  an  e-approximate  fixed  point  of  f. 

5.2.  We  construct  a piecewise  affine  homotopy  H:  S * (0,1]  — ► S - S, 
We  initiate  the  route  at  an  equilibrium  point  (xn,0)  of  the  initial 


assumption  similar 


be  ardfically  obtained,  as  Cor  example  in  [27]. 


approximations  a 
Cx,t)  e S x [0,1]  n G° 


le  with  the  fixed  points  of  piecewise  affine 
Define  H:  S * [0,1]  — ► S - S as  follows:  For 


e g(x,t)  t f(x.t)  n 


for  t 


Note  that  g(x,t)  is  an  arbitrary  vector  in  f(x,t)  n int  S 
nonempty  by  5.3.S. 

piecewise  affine  extension  of  H (on  (?)  above,  using  *.3 
5.3.7  Algorithm: 


Given  f,  N and  6 as  above,  the  algorithm  is  initiated  at  (x_,0) 
the  unique  aero  of  H in  S x [o].  Let  Tn  c S x {0}  be  the  unique 
n-simplex  of  G which  contains  Xq  (or  its  perturbed  value,  x„  + X[e] 
as  in  5.2.8)  in  its  interior,  and  0_  the  (n+l)-simplex  containing  T-. 
The  algorithmic  route,  constructed  by  a procedure  identical  to  that  in 

S x [0,1]  precludes  a "ray"  termination;  Assumption  5.3.5  and 
construction  of  H implies  that  H *(0)  n {(bd  S)  x [0,1]}  ■ $,  thus 
preventing  a boundary  termination  in  (bd  S)  x [0,1].  Thus  the 
algorithm  terminates  in  S x (1)  in  a finite  number  of  steps  when  G is 


When  the  dynamics  of  deformation  is  not  specified,  and  only  Che 

be  constructed  artifically  using  homotopy  theory.  By  using 

Remark  5.2.10  possible  equilibrium  paths  of  transition  from  the  initial 

When  f:  S * [0,»)  — ' C(S)  defines  the  evolution  of  the  economy 
an  approximate  equilibrium  path  of  evolution  4*:  (0,®)  --*■  S * [0,®) 
can  be  generated  (up  to  any  desired  value  of  t)  using  the  algorithm 
on  a refining  triangulation  of  S x [0,“)  or  S * (0,1].  For  details, 
see  Remarks  5.2.11  and  5.2.12. 

A special  case  of  the  algorithm  can  be  used  to  compute  the  unknown 
equilibrium  point  of  a given  static  economy,  lhe  details  are  similar 
to  that  in  Remark  5.2.13.  Note  that  this  special  case  reduces  to  the 
algorithm  using  homotopy  principle  of  Eaves  [30,  321,  Eaves  and  Saigal 
[34],  Merrill  [91],  etc.  Thus  our  algorithm  also  provides  a 


retrogression  of  the  algorithmic  path  c. 


By  letting  e < 


Theorem  2.4.4.  This  follows 
[0,1]  and  f(x,t)  being  nonempty 


provides  a constructive  proof  of  Theorem  2. A. 4. 


Algorithm  fi 


In  this  section  an  algorithm  ii 

with  production  (discussed  in  3. A)  under  deformation.  The  underlying 
strategy  of  the  algorithm  is  the  use  of  homotopy  technique  as  in  the 

are  quite  different  and  more  involved.  One  way  of  viewing  the 
algorithm  is  as  a parametria  version  of  Scarf's  algorithm  [107]  when 

parametria  versions  of  various  problems  to  which  Scarf's  Theorem  A. 2. 3 

parametrically  changing  Walrasian  economies. 

Recall  that  equilibria  of  a Walrasian  economy  (dfw(A)  are 
characterized  by  Definition  3 . A . 2 , as  (p ,y)  such  that 


d(p)  = 


f equilibrium  activity  levels. 


: equilibrium  prices  y 
d the  aggregate  demand 


5.A.1  Problem; 

The  problem  is  to  devise  a 
equilibrium  path  Cp,y) : [0,1]  - 


approximate 

(p(t).y(t)) 


As  in  tlie  earlier  sections  5.2  and  5.3,  a piecewise  affine 
homotopy  H will  be  defined  on  an  appropriate  subdivided  manifold.  But 
the  construction  of  the  manifold  and  the  homotopy  H are  quite  different 
from  that  of  the  earlier  algorithms,  and  in  comparison,  more  involved. 
5. A. 2 Construction  of  the  subdivided  manifold: 

The  subdivided  manifold  we  use  is  (r”  a [0,1],  (5),  where  6,  a 
subdivision  of  R?  a [0,1),  is  constructed  as  follows.  First  a 
triangulation  G of  S a [0,1)  is  constructed.  For  any  simplex  a e G, 


..,ym)  w 

. vertex  y^  " (x^,tj 
0)  is  generated;  th 
{r(yj),  r(y2),  ....  r(ym)).  The  i 


■ l(Ax 


It  is  easy  to  see  that  (R^  * [0,1],G)  is  a subdivided  manifold. 


He  start  with  H:  G^  -♦  Rn  (defined  on  the  vertices  of  G).  First 

4.3.2.  Then  for  each  point  (x,t)  e S r [0,1],  set  H(Ax,t)  ■ AH(x,t) 
all  ABO,  thus  extending  H to  r"  * [0,1]. 


H is  piecewise  affine  with  respect  to  G. 

H(5u  + (1  -5)v)  - 4H(u)  + (l-4)H(v). 

Let  c = 5a  + (l-5)b  where  u ~ ax,  v » by,  a,  b £ 0.  x,  y e 0 e G. 


(l-6)v) 


- „{c(^.xt  Ck«ib  y), 

- 

. ,H(x)  + £ii=^.H(y) 

- 6a  H(x)  + (l-6)b  H(y) 

- 6H(u)  + (1-6) Hv. 


From  construction  5.6.3  it  is  clear  that  specifying  H on  G® 
uniquely  defines  a piecewise  affine  map  H on  R,  * (0.1]  (which  is 
piecewise  affine  with  respect  to  the  subdivision  G).  In  many 
applications  of  the  alogrlthm  the  specification  of  H on  G®  will  be 


* H:  r”  * [0,1]  — ' Rn  and  a regular  value  w in  R,  the  required  route 
will  be  constructed  in  H *(w)  starting  at  in  r"  * f 0)  and 
terminating  at  a point  in  R+  * 111.  From  such  a path 
a:  [0,1]  R . * [0,1]  we  will  construct  the  required  equilibrium  path 
(P,y)t  10,11  -*  S * r£  * [0,1]  using  5.4.9. 

Together  they  guarantee  the  termination  of  the  algorithmic  path  in 

R?  * (1). 


Now  we  have  all  the  tools  required  to  set  up  the  algorithm  of 
economy  specified  as  (d(’,t),  w(t),  A( t ) ) at  time  t,  for  all  t c [0,1]. 


A boundary  assumption  is  made  similar  to  5.2.6  and  5.3.5  without 
involving  a loss  of  generality.  We  assume  Ha  = w has  no  solution  in 
{bd  r")  x [0,1],  otherwise  we  can  imbed  the  problem  and  obtain  the 
condition  artifically  by  specifying  H on  {bd  r")  x (0,1)  suitably. 

Assumption  3.4.1  we  made  for  the  Walrasian  model  is  assumed  to 


hold  for  all  t e [0,1],  i.e.,  {y:  A(t)y  + w(t)  i 0,  y i 0]  is  bounded 
for  all  t e [0,1].  This  mild  assumption  is  the  usual  one  stipulating 

nonnegative  net  supply  of  all  commodities  is  bounded. 


algorithm.  For  any  t t [0,1],  by  a suitable  intrepretation  of  a(t), 
we  will  extract  the  required  equilibrium  configuration  at  t, 

<p(t),  y(t»;  see  5.4.9. 

5.4.8  Construction  o'f  the  homotopy: 

<r;  x 10,1] , G)  is  a subdivided  manifold  constructed  as  in  5.4.2. 

G is  a triangulation  of  S x [0,11  which  generated  G.  For  any  (p,t)  eG^, 
H(p,t)  is  prescribed  as  follows. 

Set  a(p,t)  ■ max  pT^(t][  ie.,  the  maximum  profit  from  all  the 
activities  j = 1,2,  ...,k.  Choose  some  t ■ l(p,t)  such  that  pA,4)  - ir(p,  t). 
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Then  define  H(p,t)  as, 

r t)  if  x(p.t)  > 0, 

H(p,t)  - -J 

L d(p,t)  if  ll(p,t)  £ 0. 
a path  in  H *(w).  The  details  of  ci 


ay  of  extracting  (p(t),  y(t))  from  a(t)  for  each  t e [0,1]. 


the  form  (pj.t). 


(5.4.10)  H(a(t))  - I l*(t)  H(ut)  - « 

handled  by  the  algorithm  even  though  a(t)  is  used  for  its  presentation. 

Now  if  t is  in  r”  * (tJ  for  some  t e [0,1]  then  depending  on  the 
definition  of  H in  5.4.fi,  Hfu^)  in  the  equation  (5.4.10)  can  be  grouped 


follows. 


H(Uj)  - HCpj.t)  - d(pt>C)  for  i 6 A c {l,  2,  . . . , n) 
ll(ui)  - H(Pi,t)  - -A£  <t)  for  1 e r e {1,  2,  ....  n} 
tt  - f(pt,t)  as  in  5.4.8. 

Using  the  above  grouping,  (5.4.10)  can  be  written  as, 

(5.4.11)  -Z  X*(t)  A.  (t)  + Z X*(t)  d(p.,t)  - w 

r n a 1 


X1 (t)  (i  e T)  approximates 
ed;  further  that  Z X1  is 
p(t)  e T approximates  the 

e that  systems  of  the  form 


Let  0 < e i 'sCjj j"n“  ) and  A be  such  that  |d(p,t)  - d(p',t)|  i e 

A ■ {X  » (X*,  X^,  ....  Xn):  X solves  (5.4.11)}is  compact. 

Proof:  Let  e and  6 be  as  given.  Recall  t is  an  n-simplex  in 

S x (t)  such  that  T ■ conv{(p^,t),  (pjtt),  ...,  (pn,t)}  e C.  Clearly 
for  any  j,  l ( ll,  !,  . . . . n)  we  have  |p^  - Pkl  i A and 
|d(pj,t)  - d(pk,t)|  < E. 


(5.4.13)  |PJ  d(Pj,t)  - pj  d(pk.t)|  i |pj|  |d(Pj,t)  - d(pk,t) | 

< |A(Pj,t)  - d(pk,t)  | $ e. 

Now,  pj  d(pj,t)  - pjw  (by  Walras 


all  k t {1,  2,  ... 


<*■<■*»  j **<*>  •;  «•!'•>  - fj  ■ * * »*<•>  »j\'o 

i A1(t)  dtPi.t)  : P]  w 


(5.4.16) 


a & of  »o8h  «k,  k « A such 


Ik  — p(t>  - p.  lj(t)  — »I(t)  and  d(pk,t)  — d(p,t),  pk  c Tk.  In  tl.e 
(5.4.18)  -£  Ai(t)  Aj  (t)  + £ Ai(t>  d(p,t>  - w 

A*(t)  >0  1-1 n. 


Sec  5(t)  = £ X1(t);  define  y(t)  e R , the  vector  of  activity  levels, 
A 


Using  the  new  variables,  (5.4.18)  becomes, 

(5.4.19)  -A(t)  y(t)  + s(t)  d(p,t)  - w 

y(0  l 0,  C(t)  i 0. 

From  the  definition  of  H,  we  have 

(5.4.20)  if  C(t)  > 0 pTA(t)  < 0,  and 
if  y>(t)  >0  7TA(t)  l 0. 

(5.4.21)  -A(t)y(t)  = w y(t)  £ 0 
Premultiplying  (5.4.21)  by  p, 

(5.4.22)  -pTA(t)y(t)  - pTw. 

Since  FTAj(t)y^(t)  i 0 if  y^(t)  > 0 by  (5.4.20)  we  obtain  a 
contradiction  from  (5.4.22)  as  its  L.H.S.  = 0 and  R.II.S.  ■ p w > 0. 
So  C(t)  > 0 implying  pA(t)  i 0 which  is  condition  (b)  of  3.4.2  for 

Further  y^(t)  > 0 Implies  pTAj(t)  £ 0 from  (5.4.20).  Thus  we  h 

(5.4.23)  pTA(t)y(t)  = 0. 

Multiplying  (5.4.19)  'by  p we  have 

(5.4.24)  C(t)pT  d(p,t)  - pTA(t)y(t)  + pTw 

Using  (5.4.23)  and  Walras  law,  we  have  ?(t)  - 1 and  (5.4.19) 
becomes 

-A(t)y(t)  + d(p,t)  » w 

which  is  condition  (a)  of  3.4.2  for  equilibrium  at  time  t. 


limiting 


y(t)  generated 


The  next  lemma  asserts  that  an  accurate  enough  approximation  can  be 
obtained  by  selecting  a fine  enough  mesh  without  a passage  to  the  limit. 
The  lemma  gives  a set  of  bounds  on  the  degree  of  approximation  obtained 


Let  a - max  £|A*(t)|  and  a be  such  that  A(t)y  + w(t)  i 0,  y £ 0 

i i,j  3 

implleSjE^  y^  £ a for  all  t e [0,1].  Let  a be  such  that  |d(p,t)|  £ a 
for  all  (p,t)  CSX  [0,11.  Set  b = aaa.  Let  C s min  {min  w*(t)}. 

Let  S be  such  that  |d(p,t)  - d(p',t)|  i £ whenever  |p-p'|  S S for  all 
t e [0,1].  Let  6 £ min[6,  ^} . 

If  G is  a triangulation  of  mesh  6 and  r be  an  n-simplex  as  in 

in  (5. A.  11)  as  follows.  ~y(t)  e R^. 


Then  (p(t>,  y(t>)  forms 
activity  levels  such  that, 

(a)  p(t)TA,(t)  i Sa 

(b)  l>(t)TAjCt)  . -Sa 
<c)  |dl(p,t)  - wi(t)  - 


the  approximate  equilibrium  prices  and 


A*(t)y(t)  | < * + t,(f  (t/w(t,),)  + a + Seal 
p(t)  w(t)  - E 


Proof:  Given  the  choices  of  e and  6,  the  continuity  of  d ( • , t } and 


Che  details  and  refer  the  reader  to  Lemma  5.2.5  of  [107]  for  t 
underlying  arguments. 

It  has  to  be  emphasized  that  the  bounds  stated  in  5. 4. 24 
worse  than  what  occurs  in  practice. 

Now  we  sketch  briefly  the  steps  involved  in  generating  th 
a:  [0,1]  — *■  r”  * [0,1]  from  which  the  equilibrium  path 


5.4 


simplex  si 

“o  £ To  c 


Let  (p(0),y(0))  be  an  equilibrium  of  the  initial  economy, 
e n-simplex  in  G containing  p(0)  = p„,  and  aQ  be  the  (n+1)- 
thac  tq  c o_.  Generate  the  corresponding 


finite  number  of  iterations. 


termination  in  fbd  r")  a [0,11,  Since  by  construction  Cl  is  the  unique 
solution  to  Ha  = w on  R,  * [0}  the  path  is  guaranteed  to  terminate  in 

r"  « fl>. 

At  each  level  t,  (p(t),y(t))  can  be  extracted  from  oft)  using 


5. A. 28  Remark : 

homotopy  construction  as  in  remarks  5.2.10  and  5.3. 
algorithm  to  compute  possible  equilibrium  paths  of 
from  the  initial  to  Che  final  economy. 


smooth  transition 


The  algorithm  of  this  section  yields  as  a by-product  a 
continuously  refining  (homotopy)  version  of  Scarf's  algorithm  [107]  as 
Implemented  using  simpllcial  subdivision.  By  artifically  specifying  H 
fas  in  5.4.8)  on  [hd  r"1  “ [0,1]  to  be  the  Identity  map  the  analogy  .to 


parametric 


parametric  setting. 


retrogression 


APPLICATIONS 


CHAPTER 


6.1  Overview 

tn  this  chapter  we  devise  a framework  for  analyzing  various  tax 
policies  in  a dynamic  general  equilibrium  framework  using  the  algorithms 
of  Chapter  5.  The  underlying  model  is  that  of  3.6  as  modified  by  the 
inclusion  of  a very  general  tax  structure.  The  algorithms  of  5.6  would 
be  used  to  compute  the  equilibrium  paths  of  economies  under  changing 
tax  regimes. 

U.S.  [108-112.  115,  129-130).  The  chief  advantage  of  such  an  algo- 

does  not  require  any  localization  or  linearization  assumptions  and 
therefore  offers  an  appropriate  tool  for  analyzing  large  nonmarginal 
changes.  In  most  of  the  alternative  models  in  economics  based  on  a 
marginal  analysis,  parameters  such  as  demand  elasticities  are  only 
evaluated  locally  in  the  current  observable  economic  situation.  Those 
models  typically  make  many  linear  or  local  assumptions  and  are  either 

general  equilibrium  framework  used  here  not  only  permits  the  analysis 


of  large  policy  changes,  but  also  is  not  restricted  t 


such  binary 


on  the  static  models  used  in  [112,  113,  and  115],  Then  we  outline  an 

change  is  made  in  the  tax  structure  of  the  economy.  Such  paths  provide 
a dynamic  framework  for  analyzing  issues  like  the  distortionary  Impact 
of  differential  taxation  on  capital  gains  [ 1 12 ] and  comparison  of  tax 
structures  which  yield  the  same  total  revenue  [1151.  These  techniques 
are  applicable  to  evaluate  the  impact  of  policy  decisions  in  general. 


The  basic  economy  considered  here  is  the  Walrasian  model  with 


production  described  in  3.4,  characterized  by  (1)  a set  of  market 
demand  functions  d:  S — ► Rn,  (2)  the  vector  w c R,  of  the  economy's 
initial  endowments  of  all  the  commodities  l = 1,  2,  ...,  n and  (3)  a 
description  of  the  technological  production  possibilities  through  a 
listing  of  the  activity  vectors  A^,  j = 1,  2,  ...,  k.  Let  the  index 


For  convenience  first  let  us  consider  an  economy  w 


As  in  [1131  we  make 


uniform 


discriminatory 


studying  the  impact  of  distortionary  or  differential  taxation  policies. 
Thus  we  start  out  with  a general  Walrasian  model  with  an  arbitrary  set 
of  ad  valorem  commodity  taxes.  The  revenue  generated  from  the  tax 
system  is  dispersed  among  ,the  individual  consumers  each  of  whom  is 
assigned  an  arbitrary  share  of  the  total,  or  retained  by  the  government 


simply  aggregates  of 


services. 

if  3.4  the  market  demand  functions  are 
lemands  d,  each  of  which  are  derived 

the  budget  is  pw.  + 0,x).  Using  the 
e demand  functions  d are  functions  of 


1 prices  and  n 


same  arguments  of  3.2 
the  commodity  prices  j 
homogeneous  of  degree 
Walras  law.  They  can  be  denoted  ai 

S - {(p,T)  € l£+1:  Ep1  + T - 1).  I 
(p,T)  as  x.  Walras  law 
(6.2.1)  pTd(p,T)  - pTd(x)  = pTw  + x . 

gross  of  producer  output  taxes). 


the  standard  simplex 
augmented  vectors 


economy  with  producer  taxes  as  (d(‘)»  w.  A,  T) . 

6.2.2  Definition: 

Now  the  equilibria  of  (d(').  w.  A,  T)  are  defined  as  follows. 

!<  - fp.T)  e S and  y e R represent  equilibrium  sets  of  prices.  . 

equal  to  market  demand  for  each  commodity  £,  (2)  p (A-!i) ^ * 0 for  all 
j * 1,  2.  ...»  k;  with  equality  if  y^  >0.  Thus  profit  Is  maximized  at 
prices  p,  as  argued  in  3.6.2. 

At  such  an  equilibrium  we  further  have  that  the  amount  of  revenue 
distributed  t equals  the  revenue  generated  on  the  production  side  of 
the  economy  at  the  equilibrating  activity  levels  y when  prices  "p 
prevail.  To  see  this,  from  condition  (1)  of  6.2.2  we  have 
(6.2.3)  ?Td(p,r)  - pTw  + pTAy 


From  (6.2.6)  and  (6.2.3)  we  have,  T = pay,  the  total  revenue  generated 
on  the  production  side  of  the  economy. 


t.2.6  Problem: 


Given  (d(',t) , w (t>,  A(t).  T(t))  the  economy  at  time  t,  for  oil 
t e (0,11  the  problem  is  to  compute  an  equilibrium  path 
0 ■ (<5, “ (x,"y,$?):  10.11  — • S " R.  ' 10,11  of  the  changing  economy. 
If  only  the  tax  regime  is  changed  then  we  consider  the  special  case 
A(t)  - A and  w(t)  - w for  pll  t c [0,1]. 

path  using  Remark  5.4.9. 

that  H 1(w)  contains  the  algorithmic  path  a:  [0,1]  -- * R^+*  x [0,1])  has 
to  be  constructed  with  care  incorporating  the  characteristics  of  Che 
present  model.  Such  a construction  follows. 

6.2.7  Construction  of  H: 


) - A*<t)lf(t),  f 


and  t c [0,1] . As  in  5.4.8  we 
triangulation  6 of  S x [0,1]  an 
» [0,1],  For 
p‘(A(t)  - B(t)) 


■ t - 1 


e activities  of  A(t)  w 


nt  (x.t)  - <p,T,t)  c G?,  set 
um  per  unit  after  tax  profit 

,t)  - pT(A(t)  - B(t))r- 


+aJ 

tor  in  Rn. 


c path  a:  (0,11  -+  R^+1  « 


Let  S c (I.  2 k 


taxes . let  B ■ U,  2 k)  ~ 6 and  y - (1,  2 n]  ~ y. 


to  special  taxation  policy  in  s< 

tbYCO 


e equilibrium  paths  ui 


As  discussed  in  its  static  version  in  (110,  112]  the  above 

incidence  effects  and  the  efficiency  costs  associated  with  the 

consider  the  following  illustrative  example  in  which  the  existing  regime 
with  differential  taxation  across  sectors  of  income  accruing  to  capital, 
is  slowly  changed  to  a tax  regime  of  uniform  taxation. 


- {5,  6}  indexes  a heavily  taxed  sector. 

6 2 = (7,  8}  indexes  a lightly  taxed  sector. 

Y * i A ) indexes  'capital'  which  is  the  differentially  taxed 


.12)  is  Che  submatrix  o 


a particular  dynamic  tax  reform  path;  at  t a 1 
can  be  computed  using  our  algorithm  of  5.6. 


To  compare  one  equilibrium  of  an  economy  to  another  we  need  a 
criterion  function  based  on  measures  of  welfare  and  optimality.  Even 
though  no  such  generally  accepted  criteria  exists  let  us  assume  the 
existence  of  such  a measure  i)i:E  — * R where  E is  the  set  of  all  general 
equilibria.  In  studying  the  effects  of  distortionary  taxation  measures 
of  efficiency  losses  have  been  attempted  based  on  a form  of  welfare 

(51-53).  We  will  not  pursue  the  details  of  formulating  such  measures 
of  comparing  equilibria. 


further 


76 


assume  that  we  can  compare  equilibrium  paths  of 


6.3  Extension  to  Consumer  Taxes 


Extending  the  current  model  to  include  consumer  taxes  is  straight- 
forward. Such  an  extension  is  extremely  useful  allowing  consideration 
of  a wide  range  of  fiscal  phenomena  like  the  operation  of  statewide 
taxes  or  subsidies  in  a national  economy,  citywide  sales  taxes  in  a 
state  economy  and  other  issues  directly  linked  to  the  consumer. 

Each  consumer  i e I faces  a vector  of  finite  ad  valorem  consumption 

on  selected  commodi ties  for  selected  consumers  in  the  economy.  For 
convenience  we  define  a vector  such  that  • p Each 

individual  i faces  an  effective  price  vector  p + , and  pays  consumer 


Let  $(p,r)  “jEj  6j(p,T) . Then  Walras  law  for  this  case  is,  for 


(6.3.1)  pTd(p,x)  + *(p,T)  - pTw  + T. 

A general  equilbrlum  can  be  characterised  as  in  6.2.2.  At  the 


ay  + iKp.T). 


1.3.3  Construction 


fact  that  when  different  tax' systems  are  being  compared  one  often  wai 

generate  the  same  real  yield  could  be  studied  using  our  algorithms. 
Some  aspects  of  such  a model  in  the  static  framework  have  been 
considered  by  Shoven  and  Whalley  [115], 


The  underlying  model  and  the  tax  structure  incorporated  seems-  rich 

this  type  of  technique  is  that  more  than  one  distortion  or  policy  change 
can  be  simultaneously  modeled.  As  indicated  in  I 111  I a natural  exten- 


personal 


financing  social  security  out  of  increased  income  taxes  or  instituting 
or  wage  subsidy  programs  would  also  be  studied  in  variations  of  the 


6.4  Equilibrium  Paths  of  Public  Goods  Economies  under  Deformation 
Another  direct  application  of  our  algorithm  of  5.4  is  to  compute 
an  equilibrium  path  for  a continuously  changing  Public  Goods  Economy. 

For  each  period  t e [0,1]  the  algorithmic  path  would  approximate 
elements  of  a certain  subset  of  Pareto  optimal  allocations  in  the  public 


(governmental)  jurisdictions,  which  provide  public  goods  and  collect 
taxes.  All  members  (and  only  members)  of  a given  jurisdiction  consume 
the  total  of  the  public  goods  provided  by  the  jurisdiction.  Some  kind 
of  proportional  wealth  tax  is  levied  on  its  members  by  each  jurisdiction 
to  raise  sufficient  revenue  to  cover  the  cost  of  public  goods  provided. 
Private  goods  are  traded  across  jurisdictions  and  production  is  not 
jurisdiction  specific.  See  also  [97,  99]. 

Our  algorithm  of  5.4  will  be  able  to  compute  approximate 


maximization  over  private  goods  bundles  subject  to  after-tax  budget 

consistent  with  decentralized  decision  making  and  correspond  to  rather 
realistic  tax  structures. 

6.4.3  Strategy: 

Richter  (97,  99]  has  applied  Scarf's  algorithm  1107]  to  compute 
approximate  equilibria  of  public  goods  economies  in  a static  framework. 
Using  a procedure  identical  to  that  adopted  in  6.2  and  6.3  we  can  cast 
the  method  of  [99]  in  a dynamic  framework  and  implement  the  dynamic 
version  using  the  algorithm  of  5.4.  The  algorithmic  output,  intrepreted 
using  5.4.9,  will  yield  the  required  equilibrium  path. 


equilibrium  paths 


INTERNATIONAL 


equilibrium  paths  In  International  markets  under  continuous  deformation 
due  to  changes  in  tariff  structure,  import  quotas,  formation  or 
expansion  of  customs  unions,  etc.  Such  a procedure  provides  a dynamic 
framework  for  evaluating  Che  impact  of  tariff  reforms,  fiscal 
harmonization  in  ooimon  markets,  gains  and  losses  from  the  formation 


different  countries  are  treated  as  different  commodities.  The 
consequent  large  dimensionality  poses  a computational  burden;  alternate 
formulations,  as  for  example,  based  on  [54]  will  alleviate  this  problem. 
In  the  following  presentation  of  the  basic  model  we  shall  not  be 
concerned  about  such  issues.  See  Remark  7.2.5  also. 


p f R+  is  the  world  price  system.  These  prices  are  net  of  any 

which  are  produced  or  initially  located  in  country  q.  y , q c Q 

Q q 

partition  {l,  2 n)  such  that  |y  | • n , n • £jn  . 


“;f  Pq  an"  ^ y 


et  *q(p.Tq)  - gq  fq(Pq.Tq)  be  the 


A,<vV+Vp-v->Vv 

.J,  V».V  .»'i.  X\ 


7.1.4  Definition: 

Tlie  international  economy  described  above  can  be  represented  by 

endowment  of  the  world.  A is  the  technology  matrix  specifying  the 
production  possibilities,  T € Rn  ^ be  the  tariff  structure  prevailing, 
i.e.,  T denotes  the  tariffs  faced  by  the  q-th  country. 

7.1.5  Equilibrium: 

Given  the  economy  (d,w,A,T)  as  defined  above,  x = (]>,t)  e R?4^, 
t R,  represent  an  international  tariff  general  equilibrium  if 

i.e.,  Ay  + w ■=  dCp.x) . 

(ii)  profit  is  maximized  at  prices  'p,  p A^  = 0,  j ■ 1,  2,  . . . , ! 
with  equality  If  yJ  > 0 (as  in  3.4.2). 

(iii)  the  tax  revenue  received  by  each  councry  equals  that 
dispersed  by  the  government  of  that  country,  i.e.. 


7.2.1  Problem: 


Given  <d(',t),  w(t),  A(t),  T(t»  the  economy  at  time  t,  for  all 


tariff  equilibria,  i.e.,  $ ■ c (x»y»4*  ):  (0,1]  -*  S * R^  x (0,1] 

x “ (P.t)  where  p e r",  t c R^.  If  only  the  tariff  structure  Is  altered 
then  we  consider  the  special  case  A(t)  = A and  w(t)  “ w for  all 
t e 10,1]. 


5.4.  as  implemented  on  a triangulation  G of  S * [0.1)  and  a 

algorithmic  path  a:  [0.1]  r"+1  x [0.1]  we  generate  the  required 

intrepretation  of  the  algorithmic  output  are  as  in  5.4.  But  the 
piecewise  affine  homotopy,  H:  R?4^  x |0,1]  — • Rn4^  and  the  regular  value 
w e Rn+<*  (such  that  H ^(w)  contains  a the  algorithmic  path)  has  to  be 
constructed  suitably. 

7.2.2  Construction  of  H: 

Let  G be  a triangulation  of  S x [0.1).  As  in  5.4.8  we  specify  H 
on  the  vertices  (P  and  then  construct  its  piecewise  affine  extention  to 

r(x,t)  = r be  an  index  such  that  n(x,t)  ■ pTA(t),  Now  define  for 

(x.t)  t <P, 


7.2.3  Remark: 

The  construction  of  the  algorithmic  path  and  the  extraction  of  the 

algorithmic  path  Q constructed  using  Algorithm  3. A. 26  in  H ' (w) 
approximates  the  equilibrium  path  we  can  proceed  as  In  Remark  5.4.9. 

The  parallel  arguments  for  the  static  model  are  given  in  [114] . 


The  generality  of  Che  tariff  structure  T(t)  gives  us  flexibility 
for  implementing  selective  tariff  reforms  as  applied  to  a subset  of  the 
commodities  in  selected  countries.  Features  like  Import  quotas  can  be 
easily  built  into  our  framework  and  studied.  Whalley  [129]  examines  the 
process  of  harmonization  on  a common  output  tax  system  among  EEC  member 
states  within  a static  general  equilibrium  framework.  For  each  of  the 


to  changes  in  inflation  rates,  exchange  rates,  production  technologies, 

llelpman  [54]  has  improvised  Scarf's  algorithm  [107,  p.  123]  to 
compute  equilibria  in  international  markets  (in  a static  framework) 
which  alleviates  the  problem  of  dimensionality  resulting  from  the 


S the  standard  simplex  in  the  space  of  primary  factor  prices  (whose 

7.2.1.  But  the  approach  of  (541  could  be  cast  in  a dynamic  framework 
in  a straightforward  manner  as  in  the  various  applications  of  Chapters 

world  trade  equilibria. 


CHAPTER 


DYNAMIC  ANALYSIS  OF  SPATIAL  EQUILIBRIUM  MODELS 


while  the  algorithm  of  5.3  can  be  adapted  for  the  class  of  spatial 
equilibrium  models  of  8.4.  During  the  last  five  years  these  models 
have  been  studied  in  a static  framework  (for  example,  see  [4-5,  79-811) 
using  simplicial  pivoting  algorithms.  These  prior  studies  have 

In  the  discussion  below  our  main  emphasis  would  be  on  the 
adaptation  of  our  algorithms  for  a dynamic  analysis  of  the  issues 

and  the  numerous  issues  which  can  be  analyzed  using  oi 
Many  of  these  possibilities  are  discussed  in  [4-5,  79- 


8.2  Discrete  General  Equilibrium  Urban  Models 

The  salient  features  of  an  urban  model  from  residential  location 
theory  is  given  below  in  an  abstract  form.  The  emphasis  here  is  on  the 


Depending  on  the  specific  urban  issue  of  concern  (for  example  urban 
transportation  changes  14],  property  taxes  15]  and  congestion  [79])  the 
corresponding  sector  of  the  following  model  can  be  elaborated  to 
incorporate  the  necessary  detail. 


ire  Indexed 


The  land  area  of  the  prospective  city  (all  of  which  i 
put  to  residential  use)  is  subdivided  into  n sections,  i 1 

distance  d from  land  section  i.  The  residents  of  the  cil 
by  r ■ 1,  2,  R.  Every  resident  is  required  to  commul 

forth  to  the  CBD;  the  r-th  resident,  m times  a day.  Let 

Thus  p e Rn+l  represents  the  existing  price  system. 

The  welfare  of  the  residents  depend  on  consumption 
consumption  of  nonland  good  x (net  of  x,  which  is  used  os  input  into 
production  of  housing  services)  and  effective  leisure  time  T (net  of 

activities,  we  have 


(8.2.2)  x “ x + x and 

(8.2.3)  T - T + f. 

The  utility  function  of  the  r-th  resident  llr  is  described  as 

(8.2.4)  Ur  - Ur(h,x,T). 

Housing  services  are  produced  according  to  the  production  function 


(8.1 


h(x,l) 


exogenous,  i.e.,  the  residents  do  not  have  land  as  initial  endowment. 
The  initial  endowment  of  the  r-th  resident  is  exclusively  nonland  good. 


where  is  the  quantity  of  land  in  section 
commuting  expenses  applicable  in  section  i.  If 


:f  llj  is  the  optimal  value 


Let  £ and  x be  the  corresponding  demand  for  land  and  nonland  goods. 
They  constitute  components  of'his  demand  vector  d e Rn+  at  the 
prevailing  prices  p.  d = j dr  provides  the  aggregate  demand  at 
prices  p.  d:  R*’+*  ---  Rn+*  the  demand  function  is  homogeneous  of  degree 
zero  in  prices.  As  usual,  we  restrict  the  prices  to  the  standard 
simplex  S * {p:  p j - 1 , p = 0 ) . If  w c Rn+1  represents  the  vector 

of  total  initial  endowments  of  the  economy  (the  first  component  that  of 

n sections)  then  we  have  the  excess  demand  function  g:  S — * Rn+1  defined 
as  g(p)  - d(p)  - w. 


Equilibria  ci 


e characterized  a: 
d or  nonlaud  good: 


algorithm  of  5.2  can  be  applied  to  this  urban  model.  Given  that 


The  earlier  work  in  the  area  has  been  in  a comparative  static 
framework  (4-5,  79 1 . A particular  strategy  for  simplicial  pivoting 
called  SandHch  method  was  used  with  integer  labeling  1731 . Our 
algorithm  subsumes  such  methods;  using  the  equivalence  of  vector  and 
scalar  labels  in  simplicial  approximation  (see  for  example  (441)  we  can 
Incorporate  alternate  scalar  or  vector  labels  in  constructing  H in 
5.2.4  and  implement  the  algorithm  of  5.2. 


can  produce  g units  of  nonland  goods.  Thus  an  additional  constraint 
is  operative  in  maximizing  (8.2.4), 


(8.1 


Dynamic  Modeling  o 


.2  will  be  briefly  indicated.  The  strategy  is  to 
emphasise  and  elaborate  the  relevant  sector  of  the  model  incorporating 
all  the  required  detail.  We  follow  (4-5,  79]. 

8.3.1  Transportation  changes: 

To  study  effects  of  changes  in  transportation  (for  example, 
government  policies,  altering  prices  or  upgrading  facilities)  the  modal 

time,  transportation  expenditure,  etc.,  has  to  be  modeled  in  greater 
detail.  Let  the  modes  of  transportation  available  be  indexed  by 
k “ 1,  2,  ....  K.  Then  T in  (8.2.4)  is  a function  of  the  modal  choice 
k,  in  addition  to  d J and  mr  as  in  8.2.1.  Similarly  Ej  in  (8.2.6)  is 
also  a function  of  the  modal  choice  k,  d^  and  mf.  Changes  in  the 
expenses  and  quality  of  various  modes  of  transportation  Induces  changes 
in  the  equilibrium  prices,  size  of  Che  city,  welfare  levels,  population 

[41. 


8.3.2  Property  taxes: 

of  the  housing  market  is  modeled  in  detail.  Reformulate  h in  (8.2.5)  as 
a function  of  land,  structure  and  the  nonland  good,  i.e.,  h (x.f.O). 

. 0,  the  structure  can  be  measured  in  terms  of  living  area,  floor  space 


reformulated  including  the  relevant  property  taxes . (8.2.6)  is 

replaced  by, 

(8.3.3)  p0(x  + Et)  + ch(l  + p)h  = Pq  wr(l  - «) 

where  c.  is  the  cost  of  unit  housing  In  section  i (including  the  price 
of  land),  6 is  the  income  tax  rate  and  .p  the  property  tax  rate.  The 


effects  of  long  run  changes  in  property  t 


like  congestion,  effects  of 
housing  markets,  effects  of 
can  be  modeled  and  studied  in 
on  more  elaborate  versions  of 


the  model  in  8.2  (see  [79]).’ 


Deschamps  [19]  has  used  a fixed  point  algorithm  to  Che  related 
problem  of  pricing  for  congestion  in  telephone  networks.  The  problem 

as  well  as  prices  and  capacities  for  a telephone  utility  can  be 
formulated  as  that  of  computing  Kakutani  fixed  points  of  an  appropriate 
U.S.C.  mapping  (see  [ 19 J for  details).  To  compute  paths  of  approximate 
equilibria  for  such  systems  the  algorithm  of  5.3  can  be  applied.  The 

util  itics. 


8.6  Spatial  Equilibrium  Models 


spatial  equilibrium  models  discussed  here  include  as  special  cases  the 
problem  (80]. 


8.4.1  Assumptions: 

The  demand  and  supply  functions  a: 
e following  properties: 


°V  i£:  1 


e collection 


(2)  When  P(S0  then  0 c Sf (p) ; when  p£  < 0 then  S£(p)  - (0). 

(3)  There  exists  a scalar  c,  0 < c < «*  such  that  whenever 
£pf  > c we  have  (f|j  Df(p)  - Sf(p)}  n B+  ■ if 

For  motivation  for  these  assumptions  see  [81]. 

Apart  from  the  markets  discussed  the  model  has  m trading  sites 
(indexed  by  j c J)  which  are  subsidiary  institutions  located  at 
different  geographical  sites.  The  prices  operative  at  trading  sites 

role  of  trading  sites  differs  greatly  among  models. 


V»>  ■ a »’<'>■  ““‘"‘r  • « ■ ,t  4<'>- 


£ sJ(p)  j.J.ta. 
i3t  € D^(p)  j e J.  i « L. 

8.4.3  Problem: 


1 

*f(x)  - xf/xn+1  *n+1  > 1/x 


follows.  For  (x.t)  e G° , Hf(x,t)  - h£(x,t)  for  1-1,1 . 

where  h£(x.t)  e D£(«x),t)  - S£(<M>0.t).  Hn+1(x,t)  — Hf(x,t). 


H:  S » (0,1)  -*  Aff  (S  - i). 

Now  using  a algorithm  5.3.7  we  can  generate  the  required  path  of 
approximate  equilibria.  For  details  of  the  implementation,  see  5.3. 
For  parallel  arguments  in  the  static  framework  using  a vector  eanchvich 
method,  see  (81). 


equilibrium  models.  They  include  the  generalized  transportation  models 
(80),  generalized  Von  Thunen  models  of  land  use  [81J,  etc.,  as 


subclasses. 
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